Hybrid quantum-classical models as constrained quantum systems 
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Constrained Hamiltonian description of the classical limit is utilized in order to derive consistent 
dynamical equations for hybrid quantum-classical systems. Starting with a compound quantum 
system in the Hamiltonian formulation conditions for classical behavior are imposed on one of its 
subsystems and the corresponding hybrid dynamical equations are derived. The presented formalism 
suggests that the hybrid systems have properties that are not exhausted by those of quantum and 
classical systems. 
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A. Introduction 

Fundamental assumption of quantum mechanics is 
that the evolution of an isolated quantum system is given 
by the linear Schroedinger equation. On the other hand, 
all macroscopic systems usually obey nonlinear evolution 
equations of classical mechanics to an excellent approxi- 
mation. The classical and the quantum theory have de- 
veloped different formalism to successfully describe inter- 
actions between systems belonging to their respective do- 
mains. Correlations between quantum objects arc math- 
ematically captured by the direct product structure of 
the Hilbert spaces. On the other hand compound clas- 
sical systems are described on the Cartesian product of 
the component's phase spaces. Attempts to formulate 
a consistent dynamical theory of interacting quantum- 
classical, commonly called hybrid, systems are numer- 
ous as is illustrated by the following rather partial list of 
references [1-8]. Current technologies arc sufficiently de- 
veloped to enable experimental studies of the interaction 
between typically quantum and typically classical objects 
E3 ' but such experiments require detailed preliminary 
theoretical models. 

In this work the framework of the theory of Hamil- 
tonian dynamical systems is used to treat the hybrid 
quantum-classical systems and to develop a description of 
the interactions within such systems which is consistent 
with the main physically justified requirements. In fact, 
it is well known [6,11-17] that quantum mechanics can 
be formalized as a Hamiltonian dynamical system with 
the corresponding phase space and with the quantum 
observables described by functions which are quadratic 
forms of the canonical variables. More general functions 
on the quantum phase space do not have any physical 
interpretation. This formalism is used in 8] to develop 
a description of the hybrid classical-quantum systems by 
treating both, quantum and classical, formally as Hamil- 
tonian systems described in the Hamiltonian language. 
The coupling between the systems is introduced some- 
what ad hoc as if both systems were classical, just be- 
cause they are both described in the framework of the 
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Hamiltonian dynamical systems. This assumption about 
the treatment of compound systems is not trivially ob- 
vious. For example, such treatment of coupling between 
two quantum systems, both separately described in the 
Hamiltonian framework, would be incorrect. In this pa- 
per we start with the total compound quantum system in 
the geometric Hamiltonian framework. The next step is 
to consider a classical limit of one of the component sys- 
tems. To this purpose we utilize our recently developed 
theory of general quantum constraints within the Hamil- 
tonian approach [15[, and the corresponding description 
of the classical limit [l(| E3 ■ The classical behavior of 
one of the components is accomplished by constraining 
the Hamiltonian evolution so that the quantum fluctua- 
tions of the would be classical degrees of freedom remain 
minimal for all times. This in effect constrains the evo- 
lution onto a manifold which is the Cartesian product 
of the quantum phase space of the quantum subsystem 
and the manifold corresponding to the coherent states of 
the would be classical one. The evolution equations for 
the interacting hybrid systems are obtained in the macro 
limit applied on the coarse-grained subsystem and the 
constrained Hamiltonian equations. The Hamiltonian 
form of the derived evolution equations of the hybrid sys- 
tem turns out to be the same as the one postulated in Ref. 
[H and therefore satisfies a list of standard requirements 
collected and tested in @ . We also provide a discussion 
of a puzzling fact, pointed out also in Q, regarding the 
physical interpretation of functions of both the classical 
and the quantum degrees of freedom. 

B. Selective Coarse-graining and hybrid dynamics 

The framework of constrained Hamiltonian description 
for the treatment of quantum systems with nonlinear con- 
straints and its application on the problem of classical 
limit was developed and discussed in [l|| [l7| and shall 
not be repeated here. Here we apply the general theory 
in order to derive consistent dynamical equations for the 
hybrid systems. Consider a quantum system composed 
of two quantum subsystems, for convenience fancifully 
called the first and the second. The Hilbert space of the 
composite is rt = Hi ®rt2 and the quantum phase space 
of the composite is denoted M. . An arbitrary vector from 
H is denoted \ip)) and the corresponding point from M. 
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has complex canonical coordinates (ip(x),i/j*(x)) which 
are expansion coefficients in a basis of and its 

dual vector. The Poisson bracket of a two functions on 
M is given by 



{h,h}M = ~\J dx 



( Sfr Sf 2 



Sf2 Sf 2 



\5ip(x) Stp*(x) 5ip(x) Stp*(x) 



(1) 

The notation for the basis {|x))} and the integral in (1) 
should be understood symbolically and could denote a 
denumerable or finite basis and summation respectively 
The Schroedinger i.e. the Hamiltonian evolution is gener- 
ated by the function H(ip,ip*) = ({ip\Hi + H 2 + H 12 \ip)), 
where the meaning of the notation is obvious. 

For the sake of simpler presentation we shall consider 
the system such that the first of the subsystems is in 
fact given by the k— fold product of the Hcisenberg alge- 
bras, that is by the basic operators (gi, pi, .. . ,pk) 
= (q,p). We shall consider a hybrid classical-quantum 
system as a system such that the the total quantum fluc- 
tuations of the first subsystem, that is the sum of disper- 
sions of the basic obscrvables (q,p) 

k 

F(X 4 ,) = (( A <?4 + (Api)|) - mm = 0, (2) 

is preserved minimal during the evolution. This condition 
represents a nonlinear constraint on the admissible states 
of the total system. The evolution of the fully quantum 
composite system must be modified in such a way that 
the constraint is respected, and to this end we use the 
method developed in [lfl Il6| . The manifold f of the 
constraint (2) is a nonlinear symplectic submanifold of 
M. locally isomorphic with the Cartesian product Ti x 
A4 2 , where Ti is the manifold of the standard Heisenberg 
algebra minimal uncertainty coherent states (MUCS) of 
the first subsystem, denoted by |a) or \q,p), and M. 2 ~ 
T~i 2 is the quantum phase space of the second subsystem. 
Therefore, at each point \C)) of f given by 



\C)) = \a)\ua) = \q,p)\u a ), 



(3) 



there are local symplectic coordinates (q,p,uj 2 (x 2 ), 
u>2(x2)) expressed in terms of \C)) as q = ((C\q\C)), 
p = ((C\p\C)) and 0^2(^2) = (x2\(q,p\C)) ■ The vectors 
{I22)} symbolize a basis in H2 not necessarily the gener- 
alized eigen-basis of some multiplication operator as the 
notation might suggest. Notice that the requirement of 
minimal quantum fluctuations set only on the first sub- 
system automatically implies that the first subsystem is 
always in a coherent state and there is no entanglement 
between the two subsystems. Also, there can be no entan- 
glement between the degrees of freedom of the first sub- 
system. No restriction on the type of states of the second 
subsystem is set by the constraint (2), so the quantum 
subsystem can be in an entangled state. 

The fundamental assumption concerning the dynam- 
ics of the putative hybrid system is that the nonlinear 
constraint (2) is preserved during such evolution. This 



ensures that the first subsystem is minimally quantum 
(as closest as possible to classical) while the second sub- 
system is quantum in nature. Thus, our proposal for the 
dynamical equations of these coupled subsystems are the 
Hamiltonian equations given by the original Hamiltonian 
plus the additional terms that guarantee the preservation 
of the constraint (2). The resulting equations will by con- 
struction preserve the minimally quantum nature of the 
first subsystem. 

The constrained manifold T is sympl ectic and in this 
case, as was explained in detail in |Ty,|T3] , the constrained 
system is Hamiltonian with the Hamilton's function given 
by the original Hamilton's function ((ip\H\ip)) evaluated 
on the constrained manifold. Therefore, the dynamics is 
generated by the Poisson bracket on Ai and the Hamil- 
tonian 

H t = ((C(iP)\H\C(i>))) = ((i>\q,p)(q,p\H\q,p)(q,p\Tp)) 
= ({ij\H a (q,pM), (4) 

where H a {q,p) = \q,p)(q,p\ ® (q,p\H\q,p). In fact the 
constrained evolution of an arbitrary function-observable 
A(tp) = ((tp\A\tp)) on the constrained manifold is obtained 
by reducing the following equation 

A(i/>) = {A(i>),H t } M 

If, ( SA^j) SH t SH t SA&) \ 
ihj X \ 5i>{x) 5ip*(x) 5i/>(x) 5^*{x) J [ ' 

on the constrained manifold T. 

For example, before reduction on T the dynamical 
equation for q = ((V>|<?|V')) and p = are given, 

by 



1 f) M 

q = -m%H a \\n ■ 



dp 



(6a) 
(6b) 



Short computation shows that the first terms in these 
equations are in fact equal to zero on the constrained 
manifold V . In fact, for an arbitrary operator A\ acting 
only in Hi one has JT ce ]|V'))|r = 0- Therefore, 

the dynamical equations for the first system's coordinates 
and momenta are 



dH t 
dp 



P 



dq 



(7) 



Let us now compute the dynamical equations for the 
functions of the form 



^2(2:2) = {x 2 \u2W)) = (x2\(q,p\ip))- 

Starting again with the equation 
" 2[X2 > ~ ihj X \6<ip(x) SiP*(x) SiP{x) 5il>*{x), 



(8) 



(9) 
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and after somewhat lengthy calculation one obtains be- 
fore the reduction on T 



ihu 2 (x 2 ) = (x 2 \(q,p\H\q,p)\w 2 
' q dH t p dH t 
2 dq 2 dp 



Ul 2 (x 2 ) 



-(x 2 \{q,p\(p-p/2M){{i>UH* 



- -(x 2 \(qM(q-q/2M)){(mH a ]m.m 

Upon reduction on the constrained manifold T the last 
two terms are annulled and the relevant dynamical equa- 
tions can be written in the form 



ikuj 2 (x 2 ) = (x 2 \(a\H\a)\uj 2 ) 



q dH t p dH t 
2~dq~ + 2~dp~ 



u 2 (x 2 ). 
(11) 

The last term of this equation implies pure phase change 
and can be gauged away resulting with 



ihu 2 (x 2 ;i;) = (i2Ka(^)|J?|a(V0>MV0>- 



(12) 



The equation (12) has the form of a Schrocdingcr equa- 
tion for the state vector uj 2 {x 2 ;ip) = {x 2 \{q,p\ipf} € 1-L 2l 
with the Hamiltonian operator (a(ip)\H \a(ip)) acting on 
H 2 and depending on q = ((xjj\q\ip)) and p = ((ip\p\ip]). 

The dynamical equations (7) and (12) can be written 
as Hamiltonian dynamical equations in local coordinates 
on the constrained manifold T by introducing the Poisson 
bracket on T for arbitrary functions on T represented in 
the local coordinates (q-,p, oj 2 , lu 2 ) as 



{/i fo}- = V ( ^--^ - d/i 



dqi dpi 



H dpi 



± fdxj Sfl &h - 5h 6fl 
ihj \Suj 2 (x 2 ) Suj 2 (x 2 ) Suj 2 (x 2 ) Slu 2 (x 2 ) 



(13) 



Thus, the Hamiltonian form of the hybrid dynamics on 
the constrained manifold T as the phase space reeds 



q = {q,H t } T , p 

Ld 2 = {u> 2 ,H t }f, 



{p,H t } T , 

4 = W^Htjf 7 



(14) 
(15) 



where the Hamilton's function H t (q,p,uj 2 (x 2 ), uj 2 (x 2 )) in 
local coordinates on T is given by (4). 

At this point we may briefly discuss the case of gen- 
eral i.e. mixed quantum states. Such a state is given by a 
positive normalized function on M. which has quadratic 
dependence on the canonical coordinates (ip,ip*). Re- 
striction of such a function on the nonlinear submanifold 
r results with a function p(q,p, u} 2 (x 2 ),uj 2 (x 2 )) which de- 
pends quadratic on (ui 2 (x 2 ) , lo 2 (x 2 )) . On the other hand, 
p(q,p,uj 2 (x 2 ),u} 2 (x 2 )) for fixed (uj 2 (x 2 ),uj 2 (x 2 )) can be 
an arbitrary positive function of (q,p) with a unit inte- 
gral over r, since (q-,p) are not a subset of the canon- 
ical coordinates on Ai, but are physical obscrvables 
q = ((ip\q\ip^), p = ((ip\p\ip))- In fact they are a subset of 



the canonical coordinates on the nonlinear submanifold 
r. In terms of the Poisson bracket on f the dynamics 
of p{q,p,io 2 {x 2 ) 1 u) 2 (x 2 )) is given by the corresponding 
Liouville equation 

p(q,p,w 2 ,u 2 ) = {Ht(q,p,uJ 2 ,uJ 2 ), p(q,p,aj 2 ,w 2 )} T . (16) 

The constrained dynamics which preserves minimal 
value of the quantum fluctuations of one of the subsys- 
tems is only the first step. The second step is the relevant 
macro-limit so that the minimal quantum fluctuations 
still present in the corresponding coherent states can be 
neglected when compared with actual values of the dy- 
namical variables. Therefore the macro-limit should be 
applied on the equations (14) relevant for the first sub- 
system. This is illustrated in the following example. 

I. AN EXAMPLE: TWO 1/2-SPINS AND A 
CLASSICAL NONLINEAR OSCILLATOR 

Consider a system of interacting equal qubits each 
coupled to the same nonlinear oscillator. The quantum 
Hamiltonian of the total system is 



2m 



H = ea\ + ea% 

„ (17) 

where V(q) is a polynomial expression in terms of q such 
that d 2 V(q)/dq 2 \ q= o = mil 2 . The constraining and the 
macro-limit will be applied on the nonlinear oscillator 
subsystem. 

The total Hamilton function of the constrained system 
is H t = ((C{i>)\H\C(iP))), where |C(^))) = \q,p)\w). The 
complex coefficients of an arbitrary u € C 4 in the compu- 
tational basis are denoted by c\, c 2 , C3, C4 and their real 
and imaginary components are the canonical coordinates 
given by (x^yi) = v / 2(Re(ci),Im(c i )), i = 1,2,3,4. The 
expectation of the spin part of the Hamiltonian (17) in 
the vector \C)) is 

H 8 = e(yl + xj-yl-Xi)+fJ,(y2y3+yiy4 + x 2 x 3 + x 1 x 4 ). 

(18) 

he expectation in a vector \C)) of the interaction part is 

H m t = Mq{y\ +yt-yl-yl + x\ J t x\ -x\- x\)/2 

+X 2 q{yj -y 2 2 +yl- yl + x\ - x\ + xl - xj)/2, (19) 

where q = (q,p\q\q,p) is the coherent state expectation of 
the oscillator's coordinate. Finally, the \C)) expectation 
of the oscillator's Hamiltonian is 



H, 



P 

2m 



V(q) 



OO 



1 h k V^ 2k \q) 
^ 2 k kl (2mVl) k 



(20) 



where we used the explicit expression of (q,p\ V(q)\q, p) 
derived in [I7j |. In the macro- limit the term containing 
h — > is zero, leading to the Hamiltonian of the classical 
nonlinear oscillator. 

The total Hamiltonian generating the dynamics of the 
five degrees of freedom (q,p) and (xj,j/i), i = 1,2,3,4 



4 



via the Hamiltonian dynamical equations (14) and (15) 
is the sum of the three functions (18), (19) and (20). 

The dynamics of the two qubits in the form of the 
Schroedinger equation (12) is given by the Hamilton op- 
erator on %2 = C 4 which depends also on the oscillator 
coordinate q = (q, p\q\q,p) 

(q,p\H\q,p) = e&l+e&l + iin 1 x &l + \ l q(T 1 z + \ 2 q&l. (21) 

II. DISCUSSION AND SUMMARY 

In summary, we have derived from the first principles 
the Hamiltonian dynamical model corresponding to the 
hybrid quantum-classical systems that has been postu- 
lated in Q. In the derivation we have started from a 
quantum system composed of two quantum subsystems 
and then we have assumed that one of the subsystems has 
and preserves the classical properties during the interac- 
tion with the quantum subsystem. This is implemented 
by the corresponding constrained Hamiltonian dynamics. 
In this way the approach adopted in Q is justified from 
the first principles, which is our main result. 

The main properties of the hybrid dynamical equations 
in their Hamiltonian form (14) and (15) have been stud- 
ied in detail in 0] and therefore need not to be repeated 
here. However, we would like to comment on the follow- 
ing peculiar property of the hybrid Hamiltonian dynam- 
ical system already analyzed in Q. Consider the Liou- 
ville equation (16) in the case that the hybrid Hamilto- 
nian H t (q,p, u>2, ^2) an( i the density p{q,p ) u>2,^>2,) both 
depend on the same canonical pair of the classical subsys- 
tem. If the density p(q,p, W2, uj^) generates a mixed state 
on the quantum subsystem then it must be a quadratic 
function of the canonical coordinates (^(a^), w^a^)) 
corresponding to the quantum subsystem. The Hamil- 
tonian is also a quadratic function of the canonical coor- 
dinates of the quantum subsystem. However, the set of 
such quadratic functions of the quantum canonical coor- 
dinates which also depend on the classical coordinates is 
not closed under the Poisson bracket (13) on f . This is 
in sharp contrast with the purely quantum case. There- 
fore, the Hamiltonian dynamical model that corresponds 



to the hybrid system must include functions of the quan- 
tum canonical variables which do not have the physical 
interpretation of quantum observables. In fact, the hy- 
brid Hamiltonian dynamical system does not preserve the 
metrical properties of the hybrid phase space T. This is 
akin to the purely classical case where the corresponding 
dynamics preserves the symplcctic structure, i.e. the sys- 
tem is Hamiltonian, but does not preserve the metrical 
properties, which are therefore not considered as a part 
of the classical system's structure. Analogously, hybrid 
mixed states, i.e. probability densities on T, must be as- 
sumed to be of a more general form than in the purely 
quantum case. Quantum mechanical average of an ob- 
servable F in the state p: F — Tr[/3F] is reproduced with 
F = J M F(X)p(X)dX using any of the probability den- 
sities p(X) with the same first moment that is fixed by 
the requirement that the quantum expectation is equal 
to F . The fact that the quantum mixed state p deter- 
mines only an equivalence class of densities p(X), those 
with the appropriate first moment, is equivalent to the 
non-uniqueness of the expansion of the mixed state in 
terms of convex combinations of pure states and is cru- 
cially quantum property of the Hamiltonian system on 
Ai. We see that in the hybrid systems even if the ini- 
tial state p{q,p,U2,ui2) generates a quantum mixed state 
on the quantum subsystem, i.e. is quadratic in terms of 
the canonical variables of the quantum subsystem, such 
a state will evolve into a probability density of a more 
general form. This fact suggest that the truly hybrid 
systems, if existent, must be considered as conceptually 
independent class and not as such whose properties are 
exhausted by the properties of quantum and of classical 
systems. 
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